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Abstract 

Abstract: In this paper we give simplified formulas for certain polynomials which arise in some new 
Titchmarsh-Weyl m-functions for the radial part of the separated Hydrogen atom on the half line (0, oo) 
and two independent programs for generating them using the symbolic manipulator Mathematica. 



1 Introduction 



Recently Fulton [T] and Fulton and Langer [2] considered the following Sturm-Liouville problem for the 
hydrogen atom on the half line, x £ (0, oo): 



y" + + „ ]y = Xy, 0<x<oo 

X 



limiy^ y,x2 J2£+i(v4ax) = 0. 



(1.1) 



(1.2) 



In pp a fundamental system of Frobenius solutions defined at x = was introduced having the forms: 

1 



(x, A) = X 



l + Y,an{X)x^ 



n=l 



-Mg^,i{-2ix^/X), p:-- 



la 



(1.3) 
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and 



e{x,X) 



1 



2e + i 



Ki{X)(l}{x, X)lnx + X ^ + y^dn(A): 



n=l 



(1.4) 



Where an{X) and dn{X) are polynomials in A, Mg^_(_i is the Whittaker function of first kind and 



Ke{X) 



{2e + l)\{2£) 



n [4AJ^ + a' 



)2E+1 



(2£)!(2£ + r 



Then a Titchmarsh-Weyl m-function was introduced in [1] by the requirement 

9{x, A) — m£(A)0(x, A) G L2{0, oo), 

which gives as 



(1.5) 



(1.6) 



m^(A) 



aA;£(A) \ log(-2i\/A) + 'J' ( 1 



i=i l^J 

is the psi or digamma function, 



2y/xJ 

ia \ { ' (2£+l)! 



]-H2i + 27 



1 



+ 



■ 2e 2^ (-£ - 

A:!(2£ + l-A;) 



(1.7) 



where ^{z) = y^IJ 
mer symbol is defined for any complex z as (2;)^ = 
degree i defined by 



2:(z + l)...(z + A: 



Euler's constant, and the Pochham- 
- 1). Here k£{X) is the polynomial of 



(2£+l)!ll 



fl. 



For ^ = 0, we define A;o(A) = 1. The function m£{X) is an analytic function in the half planes Im A < 
and Im A > 0, has poles on the negative x-axis at the eigenvalues of the problem (jl.ip - (jl.2p . and a 
branch cut, corresponding to the continuous spectrum, on the positive real A-axis. In [2] a Pick-Nevalinna 
representation of rriQ^X) was obtained for i = and for ^ > 1 it was shown that mi{X) has a Q-function 
representation which puts it in the class N,^ of generalized Nevalinna functions with k = [^^] • 

Our purpose in this paper is to show that the last two terms in mi{X) can be decomposed into real and 
imaginary parts as 



aki{X) < 



E- 



+ 



(iVX) 



2e+i 



J 



^ (2^ + 1)! 



E 



2VA 



fc=0 



kl{2i + l-k) 



iVXki{X) + 



2^+1' 



(1.9) 



where r£(A) is a polynomial of degree £ in A. 

Effectively, the decomposition (II. 9p into real and imaginary parts becomes a defining equation for ri{X). 
Our first Mathematica program enables the decomposition to be verified. Next, in Section 3 we separate the 
left hand side of (11. 9p into real and imaginary parts by introducing polynomial representations of {—£ — t)k 
in t and 

,2 



3 



SO that the real part can be represented as a polynomial of degree £. This yields a real, somewhat explicit, 
representation for r^/(2£+l), and using it a second Mathematica program shows that this real representation 
yields the same result as the first Mathematica program. 



2 The First Mathematica program 



A simple program in Mathematica can be used to verify that the polynomial r£(A) in p.9|) is real valued. 
This makes use of the built-in function Pochhammer [a, k] which executes the multiplications in the 
Pochhammer symbol (a)fc. 

Using the second expression in (jl.Sp for K£{X), we have using (jl.Sp that 



-aki{X) 



Ke{X) _ {-2iVXf'+H-i-Phe 



+1 



2i + l 



[{2i + r 



112 



Accordingly, solving ()1.9p for 

n{X) . 1 



21 + 1 



E 



J 




where c : = 



21+1 



(2^ + 1)! 



■ 21 2^ (-£ - 

k\{2l + \-k) 



(2.1) 



Following is the Mathematica program which implements equation (12. ip . The output for "^^t+i 
i = 1,2,3, and 4 is shown. The program was executed up to ^ = 30 showing that ri{X) remained real- 
valued. Observe that the constant a multiplies all terms of which in the Mathematica output is ansl. 
This is also poved in (j3.7p below. 



First Program 

Program input: i 

£ = 4; 

While[£<6,/3=^; 

b = Expand[(-2/^)^'+i Pochhammer [-£- ft 2£+i] . ^^^^.^ generates Ke/{2e + 1). See (1.5).*) 

„ _ AnartfSimnlifvfi^^^^ Pochhammcr[-^-/3,fc]2fe ni 

c — Apart [simpmy[ (^2e+iy. 2^k=o fci(2f+i-fc) JJ' 

Print ["For £ =", £]; 

ansl = Simplify[c + e + /]; 

Print["ansl = "ansl]; (*This generates the RHS of (2.1).*) 
ans2 = Series[ansl, {A, 0, Print[" = ",ans2]; 

ans3 = Simplify[5 + c + e + /]; Print ["ans3 = " , ans3]; (*This generates the RHS of I^M*) 
ansA = Series[ans3, {A, 0, £}]; Print [" = ", ans4]; 
£ = £ + 1] 

Program output: 

For ^ = 1 

aA 

~36 

— (Sa^ + lOaA) = \ 

72^ ^ 24 36 

(2.2) 



ansl = 
ans3 = 



A 



For £ = 2 



ansl 



aA(a2 + 13A) a^X ISaA^ 



7200 7200 7200 



25a^ + 476a3A + 1288aA2 llda^X leiaA^ 

ans3 = = 1 1 

172800 6912 43200 21600 



For ^ = 3 



ansl 
ans3 



aA(a^ + 46a2A + 400A2) _ a^A 23a3A2 aX^ 

8467200 ~ ~ 8467200 ~ 4233600 ~ 21168 

490^^ + 26840^ A + 35656a3A2 + 88896a A^ 

508032000 

a"^ 671a3A 4457a3A2 463aA3 



For e 



ansl 



ans3 



10368000 127008000 63504000 2646000 



aXja^ + 107a^A + 3124a^A^ + 22548A-') 
32920473600 

a^A 107a5A2 781a3A3 1879aA4 



32920473600 32920473600 8230118400 2743372800 
761a^ + 90200a^A + 3204208a^A^ + 36438400a^A^ + 86960256aA'^ 

36870930432000 
761a^ 451a^A 4087a3A3 226459aA'' 



36870930432000 184354652160 11522165760 96018048000 



Output for the polynomials 



which arise in the representation of the nii function from [H Equation(8.15)], 

me{X) = ke{X) |-alog(-2f\/A) - a^- ^1 - ■^^=^ - 27a + ix/Xj + Pi{X) 
is also given below for ^ = 1,2, 3, 4. 



Pi = TTT + 



a'^ 5aA 



24 36 
a^ 119a^ 161a 2 
~ 6912 ^ 43200 ^ 21600 ' 



a 



7 «71r,5 



671a5 . 4457a3 ,0 463a 



10368000 127008000 ^ ^ 63504000 ^ ^ 2646000 ^ 

761a9 4510"^ , 4087a^ ,2 113870^ ,0 226459a 

Pi = o^o>.nnon.onnnn + 7177^7777:^777777 A + , ^^^^^ , ^^^^ A^ + 777777777777777 ^ + 



36870930432000 184354652160 47029248000 11522165760 96018048000 
Here pi is printed as ans3 in the above program. 



3 Explicit Representation for 

In this section we give a method for separating the expression in equation p.9|) into real and imaginary 
parts, yielding a real representation for the polynomial r^(A). The difficulty arises from the complicated 
product in the Pochhammer symbol {i — P)k where f3 = Replacing /3 by a real variable t, we let the 
coefficients of the polynomial {—£ — t)k be defined by 

fe— 1 k fci k2 

9kit) := i-i-t)k = Hi-i + j-t) = Y,aik,n)e = ^ a(A:, 2j)t2j' + ^ a(A;, 2j + (3.I) 

j=0 n=0 j=0 j=0 

where ki = [|] and ^2 = [^y^]- Here a{k,n) = a£{k,n), and we are interested for fixed £ to have a£{k,n) 
available for all < A; < 2£, and all < n < A;. For n = 0, the constant term is ai{k, 0) = nj=o + j) ~ 
{—i)k, k = 0,1, ■ ■ ■ 2i. Formulas for larger values of n become increasingly more complicated and are not 
known in closed form. We can, however, represent the real and the imaginary parts of (jl.Op in terms of 
a{k,n). Putting t = P = (f3l^ gives 



n=0 



g,{P) = Y.a{k,n)(^yi-{X--2r 



j=0 



:X;(-i)Mfc,2, + i)(|)''^'(A-- 



2j+l 



(3.2) 



j=0 



Now for the (i\/A)^^"'"^ term in (|1.9p we have (i)^^"^^ = (— i = 0, 1,2, ... Thus the second sum in 



!jl.9p may be written as 

21 



k=0 

2e 



{2e + i)\ki{2e + 1 - k) 



5fc(/5)(A 



^)2e+i 



^ {2£ + l)\k\{2£ + 1 - k) 
2« ofc 



X;(-iya(fc,2j + l)(|)'^'^\A^r+i-(^^+i) 



j=0 
ki 



k=0 



(2£ + 1)!A:!(2£ + 1 - k) 



2M, 



2£ 



j=0 

2^Bk 



21 

^ (2Z+ 1)!A;!(2^ + 1 - A:) ' " ^ {2£ + l)\k\{2£ + l-k) 

k — A? — 



(3.3) 



where 



and 



k2 



a\2i 



^fc = (-i)'+^E(-i)Mfc'2j + i)U) 

j=0 

B, = {-lYj:i-iyaik,2j){^y\'-^x'^ 

j=0 



(3.4) 



(3.5) 



Similarly, we may decompose the first term in (|1.9p as 



E 



1 



m 



E 



4Am 



^ ^ 4Am2 + a2 

2VA m=l 



+'E 



2aA2 



(3.6) 



»Tt=l 2\/A m=l m=l 

Using (|1.8|) and cancelling one 4Aj^ + factor for each term in the above sums, we find that 



^ 1 

y — — 

.m=l 2VX. 



[(2£+l)!]2 



m=l j=l 



m + 1; 



112 



m=l 



(3.7) 



Combining (j3.3p and (j3.7p and taking real and imaginary parts of (jl.Op . we thus obtain: 



\/AA:£(A) 



^aVA_ 



m=l j=l 



2^ 



+E 

fc=0 



(2^ + l)!A:!(2£ + l-fc) 



(3. 



and 



4Aa 



n{\) 

2i + l [(2^ + l)!]2 



m=l j=l 



2e 



+E 



^ (2£ + l)!A;!(2^ + l- A:)' 



(3.9) 



Since every involves a as a factor, it follows that a factor of a occurs in both terms on the right hand 
side of (|3.9p . Hence a = will give ri{X) = 0. 

To bring the expression in (13.91) into a simpler form, we now isolate the coefficients of the powers of A. 
For the second term in (j3.9p we insert the formula (|3.4p for Ak into the sum, change the summation index 
(m = i — j), and then interchange the order of summation: 



2''Ak 



2i 

^ (2Z+ 1)!A:!(2£ + 1 - k) 
^ ^ 2'=(-l)^+i+^a(A:, 2j + 1) (f )'^+' A^"^' 

A:=l j=0 



{2e + iy.kl{2£ + 1 - k) 



21 i 

E E 

k=l m,=£—k2 

e 21 

E E 

m=l fc=2(f-m)+l 



a\2{i-m)+l 



2fc(_l)™+ia(A:, 2(^ - m) + 1) (§) 
{2l + l)\k\{2l + l-k) 

2*^(-l)™+ia(A:, 2(^ - m) + 1) (|)2(^-'»)+i 
{2£ + l)\k\{2£ + I - k) 



(3.10) 



where 



di := - 
^ 2 



2£ 

E 

fc=2(f-i)+l 



2fc(_iym+i (i^A* 1) a{k,2{i-j) + l) 
(21 + iy.k\{2i + 1 - k) 



(3.11) 



Similarly, to isolate the powers of A in the first term in (j3.7p we first define the coefficients 7(m, n) of the 
i — 1 degree polynomial, 

n^ + ^)=E7Kn)A". (3.12) 

j=l ^ ^ n=0 

Insertion of this into the first term in (j3.9p . changing the summation index (j = n + 1), and interchanging 
the order of summation then yields: 



-4Aa 



WTTW 

[(2^ + l)!]2 
[(2^ + l)!P 



^mn(4Aj2 + 

m=l j=l 



EE 

.m=l ?i=0 

EE 

m=l j=l 

£ e 

EE 

j=l m=l 



7(m, j - 1) 



m 



- 1) 



m 



A^ 



A^ 



(3.13) 



where 



a4^(^!) 



m=l 



7(m,i - 1) 



m 



(3.14) 



Putting p.lOp and ()3.13p in p.9p we thus have for the polynomial r£{X)/{2i + 1) the representation. 



2£ + 



(3.15) 



Here it is clear that a is a common factor in Cj and dj for all j, and hence for the case o = 0, = 0. The 
expression (|3.15p is somewhat more explicit than the expression (j2.ip and (j3.9p since the coefficients of the 
powers of A are isolated, and there are no complex terms present. On the other hand, further simplifications 
are certainly desirable; however, this requires closed form formulas for a{k^ n) and ^{m, n) which remain 
elusive. However, the formulas for Cj and dj are easily implemented using a symbolic manipulator. 



4 The Second Mathematica Program 

As an independent check on the first Mathematica program we implemented the formulas (j3.1ip and (|3.14p 
for dj and Cj and computed the polynomial of degree £ in ()3.15p . This required computing and storing the 



coefficient a{k,n) and ^{m,n) in (j3.ip and (j3.12p . Following is the Mathematica program that does this 
to compute r^{\)/{2l + 1). The output for I = 1,2,3 and 4 is shown. The program was executed up to 
£ = 30 and gave exact agreement with the first Mathematica program. 



Program input: I 



Second Program 



^ = 4;A; = 0; 
While [A: < 2£, 

p\tJ\ = Expand[Pochhammer[— ^ — t, fc]]; n = /c; i = 0; 

While [i < n, 
a[k,i\ =p[0]; 

p[t_\ = Simplify[Expand[(p[t] - a[k,i]) /t]\; 

k + +] 
Clear [A:] ;j = 1; 
While [j < 

"UJ — 2^fc=2(^-i)+l {2l+l)\k[{2l+l-k) iJ-r-^J 

Clear [j] ; m = 1; 
While [m < 

<i[t-] = {WU + ^) ) / + 4^) ; - = ^ - 1; ^ = 0; 

While [i < n, 
7[m,i] = g[0]; 

q[t.] = Simplify[Expand[(g[t] - 7[m,i])/t]]; 

i + +] 
m + +] 
Clear = 1; 
While [j < 



^Uj — ((2^+1)!)^ 2^m=l m + 

PrmtE5=i(c[i]+t^[i])t^] 



Output of second program 

For £ = 1 : 



For 



at 



{2i + 1) 36 

ah Uat^ 



For 



For I 



{21 + 1) 7200 7200 

a^t 23ah'^ at^ 



(2£ + l) 8467200 4233600 21168 
re a'^t WTaH'^ 781ah^ 1879ot^ 



(2^ + 1) 32920473600 32920473600 8230118400 2743372800 

This output is in agreement with the output of the first program, and was also checked up to ^ = 30. The 
two Mathematica programs verify that the imaginary part of the expression (|2.1|) for ri{X)/{2i + 1) is zero 



n 



for all values of H. for which the programs were executed. A general proof that this is true for all that 
is, a proof of (|3.8|) (or, equivalently, a proof of (|c{.9|) ) requires that the a{k,n) and 'y{m,n) coefficients be 
obtained in a simplified form, and this appears to be a formidable task. It is quite difficult to perform 
induction on £ to prove (|3.8p or (j3.9|) because of the complicated manner in which a{k,n) = ai{k,n) and 
7(m, n) = j£{m, n) change with £. Nevertheless, it may be possible to construct rigorous proofs by making 
use of some combinatorial analysis. For example there are general forms for the coefficients of a polynomial 
having k known roots Xj, i = 1,2, ...,k. Sen and Krishnamurthy [3], for example, show that 

k k 

j=l m=l 

with 

am ■■= i-irSm ■■= (-ir Yl ^l^2..-Xm, (4.2) 

where the sum is taken over all the products of Xi taken m at a time. Thus, for the polynomial gk{t) of 
(|3.1|) we have 

k-l k / \ 

9k{t) = {-!)' Hit - {j - £)] = i-l)'ll[t - (j - 1 - i)] = {-!)' (t'+Y. ^rnt'-"^ ' (4-3) 

j=0 j=l \ m=l / 

with 

a„ = a£(/c,m) = (-l)™^tit2...tm (4.4) 

where the sum is taken over all the products of tj := j — 1 — £ taken m at a time. But, a general solution 
for ttm as a function of k and £, valid for all £ and all < A: < 2£, remains elusive. Another idea would be 
to establish (j3.8p or (j3.9p by induction on £, but it unfortunately appears difficult to employ the induction 
hypothesis. So, an analytic proof of (|3.8p or (|3.9p valid for all £ > 1 remains an open problem. 

5 Conclusion 

The simplification of a^(A;, n) and a rigorous proof of (j3.8p and/or (j3.9p remain as open problems. We have, 
however, given in this note two Mathematica programs which implicitly establish both of these results. 
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